A robust all-non-optical method for the characterization of single-shot few-cycle laser 

pulses 



Zhangjin Chen, 1 T. Wittmann, 2 and C. D. Lin 1 

1 J. R. Macdonald Laboratory, Physics Department, 
Kansas State University, Manhattan, Kansas 66506-2604, USA 
2 Max-Planck-Institue fur Quantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany 

(Dated: June 12, 2009) 

An all non-optical method for accurately determining the pulse parameters of individual few-cycle 
laser shots is presented. By analyzing the "left" and "right" asymmetry of high-energy photoelec- 
trons along the polarization axis using the recently developed quantitative rescattering theory, we 
show that the carrier-envelope phase, the pulse duration and the peak intensity of each single-shot 
pulse can be readily retrieved. Unlike optical measurements, this method provides the laser intensity 
in the interaction region directly. 
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Observation of the interaction of laser pulses with mat- 
ter requires the knowledge of the waveform of the laser 
field. The waveform is characterized by its pulse enve- 
lope, including peak intensity and pulse duration, as well 
as the carrier-envelope phase (CEP) which measures the 
offset between the peak of the electric field and the peak 
of the envelope. Specifically such a CEP-fixed waveform 
can be written as E(t) = Eo(i) cos(ujt + </>), where ui is 
the frequency of the carrier wave, <j> is the CEP. 

The stabilization of the CEP is of great importance, for 
example, in the generation of isolated attosecond XUV 
pulses [l|, and in controlling the dissociation dynamics 
of molecules For atomic targets, the CEP affects the 
relative yield of high-energy above-threshold-ionization 
(HATI) electrons along the left vs the right directions 
of the polarization axis, and the ratio of the left/right 
photoelectron yield has been used for the determination 
of the CEP of few-cycle pulses In such measure- 
ments the CEP of a laser beam must be stabilized. How- 
ever, the stabilization of the CEP is rather complex. It 
is often determined by averaging over the jitters from 
many laser shots, thus limiting the temporal resolution 
of "CEP-stabilized" pulses. Recently, Wittmann et al. Q 
reported that the left/right photoelectrons from individ- 
ual shots can be measured. In this Letter, we show that 
based on the recently developed quantitative rescattering 
(QRS) theory [|, the CEP, the pulse duration, together 
with the peak intensity at the laser focus, can all be ac- 
curately determined simultaneously from the measured 
electron momentum spectra, thus providing a complete 
non-optical method for the full characterization of each 
few-cycle laser shot. The method is robust and the laser 
pulses can be characterized within a few hours using typ- 
ical personal computers. 

From the QRS [f|, the momentum distributions of 
HATI electrons can be expressed as 

D(p,6) = W(pr)a( Pr ,6 r ) (1) 

where W{p r ) is the returning electron wave packet and 



a(p r , r ) is the elastic differential cross section (DCS) be- 
tween the field-free electrons and the target ion. The 
QRS is a quantitative version of the well-known rescat- 
tering model [1, Q ■ In this model, electrons which are 
released earlier in the laser field may be driven back to 
recollide with the target ion. In Eq. ([T]), the momentum 
p r of the returning wave packet (not directly measured) 
is related to the measured photoelectron momentum p 



pcosO = —A r TPrCos9 r , 
p sin — p r sin r 



(2) 
(3) 



where A r is the vector potential of the laser field at the 
time of recollision. In Eq. (J2]), the upper sign refers to the 
right side while the lower sign refers to the left side. (All 
the equations are written in atomic units unless other- 
wise noted.) Following the QRS [IS], the magnitude of 
A r =p r /1.26. The angles 9 and 9 r are defined with respect 
to the laser polarization axis, and 9 r is measured from 
the direction of the "incident beam" , i.e., of the returning 
electrons. For HATI electrons along the polarization axis, 
0=0 (right) or tt (left), 6> r =7r, and p=2.26p r /l. 26=1. 79p r . 

According to the QRS @, i], the wave packet W(p r ) 
depends mostly on the lasers only, while the DCS is 
solely the property of the target. Due to this separabil- 
ity the experimental HATI spectra can also be written as 
the product of a volume-integrated "wave packet" W(p r ) 
with the DCS [B| . For each pulse, there are two returning 
wave packets, one from the left, the other from the right, 
along the polarization axis. Fig. 1(a) shows the "left" and 
"right" electron spectra from a typical single-shot mea- 
surement. Using Wn(pr) = D(p,9 = 0)/a(p r ,9 r = it), 
where a{p r ,9 r = 7r) for atomic targets are easily calcu- 
lated, for example, the right wave packet can be obtained. 
Similarly a left wave packet W l (Pr) can also be obtained. 
By comparing these "experimental" wave packets with 
those W(p r ) obtained theoretically the laser parameters 
used in the experiment are retrieved. Since the wave 
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packet in the QRS is obtained from strong field approx- 
imation, the calculation is a few thousands times faster 
than from solving the time-dependent Schrodinger equa- 
tion. This speed-up makes it possible to carry out near 
real-time retrieval of laser parameters from experimental 
data. To perform volume integration, we assume that the 
spatial distribution of each laser shot is Gaussian and all 
the electrons from the interaction volume are collected. 
Thus only the peak intensity at the laser focus and its 
pulse length (defined as the full-width at half-maximum 
of the intensity) , in addition to the CEP, will be retrieved. 
We also assume that the peak intensity and pulse dura- 
tion of each individual shot do not change, i.e., only the 
CEP is changed randomly. The wavelength is assumed 
to be known. 

In Ref. 4500 single-shot data were collected. We 
first determine the peak intensity and pulse duration used 
in the experiment. For this purpose, we define a single 
quantity called energy moment M for each shot, 

CW(p r )dp r W 

where the momentum p r of the returning electron is re- 
lated to the momentum p along the polarization axis (and 
energy Ei = pf/2) by p=1.79p r . For example, from the 
experimental electron spectra like Fig. 1 (a) (or better the 
spectra summed over all the shots, on the left or on the 
right) we estimate E\ and E 3 , where E\ is close to about 
5U P and E 3 is about 10U P , where U p is the ponderomo- 
tive energy. These selections are made since the QRS is 
valid only for HATI electrons. The precise values of E\ 
and £3 are not important. Using the left wave packet for 
all the 4500 shots, 4500 values of M's are calculated from 
Eq. ([!]), using fixed E\ and E 3 . These calculated values 
are displayed in Fig. 1 (b) where the horizontal axis is di- 
vided into 90 sections. The moments M calculated from 
the first 50 shots are placed in the first bin, at the verti- 
cal positions corresponding to the values of M. The M's 
from the next 50 shots are placed in the second bin. The 
process continues till the energy moments from all the 
shots are registered. Note that the M's are distributed 
nearly uniformly within a band. The average value of M, 
or M, was calculated to be 16.46 eV and was found to be 
independent of the pulse duration. Using the QRS, the 
M values (averaged over the whole 2ir range of the CEP) 
were found to be 15.64, 16.05 and 16.46 eV, respectively, 
for peak intensities of 1.2, 1.3 and 1.4 x 10 14 W/cm 2 . 
By choosing peak intensity at 1.4 x 10 14 W/cm 2 , and 
for pulse durations of 4.5, 4.7 and 5.0 fs, respectively, we 
found that the best fit to the (vertical) band width is for 
pulse duration of about 4.6 fs, see Fig. 1(c). Thus the en- 
ergy moments calculated over the single shots allow the 
determination of the peak intensity and pulse duration 
easily. 

In the above analysis, only the energy moments from 
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FIG. 1: (Color online) (a) Typical single-shot left- (solid line) 
and right-side (dashed line) electron energy spectra along the 
polarization axis, (b) Energy moment of left-side spectra from 
experimental measurements (dots) compared with theoretical 
calculations at peak intensities of 1.2 (dash-dot line), 1.3 (bro- 
ken line), and 1.4 x 10 14 W/cm 2 (solid line) with pulse dura- 
tion of 4.5 fs; (c) Same as (b) but for theoretical calculations 
at peak intensity of 1.4 x 10 14 W/cm 2 , and pulse durations 
of 5.0 (dash-dot line), 4.7 (broken line) and 4.5 fs (solid line), 
respectively; (d) Energy moment of right-side spectra from 
experimental measurements. The experimental data are from 
Wittmann et al. 4] measured for single ionization of Xe in 
laser pulses with the wavelength of 760 nm. In (b-d), experi- 
mental data from all 4500 shots are shown, see text. 

the left wave packets were considered. If the right and 
left detectors are exactly identical, then the same peak 
intensity and pulse duration should be obtained from the 
right wave packets. From Figs. l(c,d), it is clear that the 
two detectors are not exactly the same. If we were to 
use the data from Fig. 1(d), we would obtain a peak 
intensity of 1.33 x 10 14 W/cm 2 and pulse duration of 
4.8 fs. We checked that these conclusions are not changed 
much when the values of E\ and E3 are varied. 

Once the peak intensity and pulse duration are known, 
we retrieve the CEP for each shot following the proce- 
dure of Wittmann et al. [4]. In their method, between 
Ei and E 3l another intermediate energy E2 was chosen 
(see Fig. 1(a)). The total electron yield Yl between E\ 
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FIG. 2: (Color online) Comparison of the asymmetry ellipse 
from experiment with theory. The experimental points are 
calculated from the data of Wittmann et al. 3. The peak in- 
tensity used in the theoretical simulations is 1.4 x 10 14 W/cm 2 . 
(a) The energy range used is (Ei, E 2 , E 3 )=(37.9, 57.5, 64.8) 
eV, and pulse durations used in theoretical simulations are 
4.5 (broken line), 4.7 (solid line) and 5.0 fs (dash-dot line), 
respectively; (b) The energy range is (Ei, E2, I*y = (38.6, 48.7, 
82.0) eV and the pulse duration for theory is 4.7 fs (solid line). 



and E 2 from the left detector is evaluated, and a Yr from 
the right detector in the same energy range is calculated. 
Define the asymmetry A x — (Yl~Y r )/ (Y L +Y R ). A sim- 
ilar asymmetry parameter A 2 is defined for the electron 
yields between E 2 and E 3 . Using (E 1 , E 2 , £ 3 ) = (37.9, 
57.5, 64.8) eV as in Wittmann et al. 4], the (Ai,A 2 ) for 
each laser shot is plotted as a point in 2D, and the results 
for all the shots are shown in Fig. 2(a). On top of the 
plot, three theoretical curves are shown, for peak inten- 
sity of 1.4 x 10 14 W/cm 2 and pulse durations of 4.5, 4.7 
and 5.0 fs, respectively. From the three curves, a dura- 
tion of 4.7 fs gives the best overall fit to the experimental 
data. This number happens to be the average of 4.6 fs 
and 4.8 fs derived from Fig. 1. Note that the theory curve 
is simply a Lissajous parameter plot of Ai and A 2l ver- 
sus the implicit variable, the CEP. The theory expects a 
perfect ellipse. Due to the intrinsic errors in experimen- 
tal electron spectra, the experimental "ellipse" acquires 
a width, and the ellipse is distorted due to the difference 
in the left and right detectors. 

The shape of the ellipse depends on the choice of 
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FIG. 3: Absolute CEP extracted for experimental measure- 
ments of Wittmann et al. [J from shot number 1000 to 1050. 



(Ei, E 2 , E 3 ) used in the calculation of and A 2 . In 
Fig. 2(b), we show another choice of these parame- 
ters. The size and the orientation of the ellipse are 
changed. However, the actual retrieved CEP's are in- 
sensitive to such choices. Take laser shots, #829, #1138 
and #4000 as examples. We retrieve the CEP by drawing 
a straightline from (A\, A 2 )=(Q, 0) to the experimental 
point (marked by large yellow dots). From the intercept 
of this line with the theoretical curve, we read out the 
CEP. For these three shots we found their CEP values 
arc 129°, 279° and 16°, respectively, with an error of 
about 3°. No effort was made to optimize the choices of 
the three energy points. 

To illustrate that the CEP of each laser shot indeed 
varies randomly, we show the retrieved CEP for shot 
numbers from 1000 to 1050 in Fig. 3. Clearly the CEP 
varies randomly from shot to shot. 

The present method can also be used to determine the 
CEP of the "phase-stabilized" laser pulses. We have elec- 
tron momentum spectra from Kling et al. [l(]| . From their 
HATI electron momentum spectra, we deduced that the 
peak intensity is 3.0 x 10 13 W/cm 2 and the pulse dura- 
tion is 5.0 fs. The mean wavelength of the laser used was 
738 nm. In Fig. 4(a), we show the parametric plots using 
(Ei,E 2 ,E 3 )= (10.0, 13.0, 20.0) eV for the 19 measured 
points. The "size" of the ellipse is about the same as in 
Fig. 2 since the pulse duration of 5.0 fs is close to the 
4.7 fs in Fig. 2. However, the scattering of the experi- 
mental data points from the theoretical ellipse is much 
larger for these "phase-stabilized" laser pulses. This large 
scattering reflects the lack of good phase stabilization. 
As shown in Wittmann et al, performing single-shot 
measurements on these "phase-stabilized" pulses shows 
shot-to-shot CEP variations up to about 20°. Using the 
present method to retrieve the CEP, we found that the 
retrieved CEP depends more sensitively on the values of 
(E\,E 2 ,E 3 ) used in the analysis, see Fig. 4(b), where 
two sets of energy values were used to obtain the CEP 
for each measurement. The straightline in Fig. 4(b) was 
drawn so that the line best fits the deduced CEP's from 
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FIG. 4: (Color online) Retrieval of CEP for "sequential" 
"phase-stabilized" measurements, (a) Comparison of asym- 
metries from experiment (crosses) with asymmetry ellipse 
from theory (solid line). The asymmetries are calculated us- 
ing (E 1 ,E 2 ,E 3 )^(W.0, 13.0, 20.0) eV. Experimental data for 
19 shots are from Kling et al. [l(j using wavelength of 738 nm. 
The peak intensity in the theoretical simulation is 3.0 x 10 13 
W/cm 2 and the pulse duration is 5.0 fs. (b) The absolute 
CEP's for some shots extracted from (a) (circles), compared 
to those extracted using (Ei, E 2 , £ 3 )=(10.0, 15.0, 20.0) eV 
(crosses). (c,d) Similar comparison for "phase-stabilized" long 
pulse measurements from Kling et al. [llfl . The data consist 
of 40 sequential shots with wavelength of 760 nm. In (c), 
shot numbers from 21 to 31 are marked. The peak intensity 
used in the theory simulation is 1.05 x 10 14 W/cm 2 and the 
pulse duration is 7.0 fs. The energies used for the asymmetry 
calculations are (Et, E 2 , E 3 )=(40.0, 50.0, 60.0) eV and the 
retrieved phases are shown as circles in (d) . For the crosses in 
(d) the energies used are (40.0, 45.0, 55.0) eV. The straight- 
lines in (b) and (d) are best fits to the deduced CEP's. 



tensity is 1.05 x 10 14 W/cm 2 , close to the values from 
Micheau et al. T2|. For these longer pulses, the asymme- 
tries Ai and Ai are much smaller, see Fig. 4(c). Due to 
the remaining jittering of the phase stabilized pulses from 
shot to shot, the scattering of the experimental asymme- 
tries around the theory ellipse are much larger. For the 
ten measurements over the 27r range, we notice again 
that the retrieved CEP's depend more sensitively on the 
(Ei, E 2 , E 3 ) used for the retrieval, see Fig. 4(d), similarly 
to what was observed in Fig. 4(b). 

In summary, using the recently developed quantita- 
tive rescattering theory, we propose an all-non-optical 
method to efficiently retrieve accurate peak laser inten- 
sity, pulse duration and the CEP of each single laser 
shot from the measured high-energy ATI electron spec- 
tra. The method is very robust. The computational 
effort is very small and requires no iterations. Unlike 
Wittmann et al. we show that there is no need to 
retrieve the peak intensity and pulse duration using the 
optical method. In the future, CEP dependence mea- 
surements can be carried out using non-phase-stabilizcd 
lasers. By simple CEP tagging, the strong-field effects on 
the waveform of ultrashort pulses can be accurately car- 
ried out to achieve temporal resolution of a few attosec- 
onds. The present all-non-optical method also character- 
izes the laser pulses in the interaction region directly and 
avoid errors introduced by the propagation of the laser 
beam. 
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successive measurements. On the whole the phase de- 
creases as the "shot number" increases, thus the phase 
has been stabilized at least partially. 

The determination of CEP is easier for shorter pulses. 
Previously, Kling et al. [ll[ reported HATI electron mo- 
mentum spectra for lon ger pulses. Their data were ana- 
lyzed by Micheau et al. [T2J earlier, also using the QRS, 
but the method assumed that CEP differences between 
successive measurements are constant. For this same set 
of data, Micheau et al. found a pulse length of 6.7 fs and 
peak intensity of 1.05 x 10 14 W/cm 2 . Using the present 
method we analyzed the 40 "shots" from Kling et al. ll| , 



and found that the pulse length is 7.0 fs and the peak in- 
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